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^ 1 We present new measurements of the power spectra of the E-mode of CMB 

polarization, the temperature T, the cross-correlation of E and T, and upper 
limits on the B-mode from 2.5 years of dedicated Cosmic Background Imager 
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ABSTRACT 



(CBI) observations. Both raw maps and optimal signal images in the w-plane 
and the sky plane show strong detections of the E-mode (11.7a for the EE power 



spectrum overall) and no detection of the B-mode. The power spectra are used to 
constrain parameters of the flat tilted adiabatic ACDM models: those determined 
from EE and TE bandpowers agree with those from TT, a powerful consistency 
check. There is little tolerance for shifting polarization peaks from the TT- 
forecast locations, as measured by the angular sound crossing scale 9 = 100/ £ s = 
1.03 ± 0.02 from EE and TE cf. 1.044 ± 0.005 with the TT data included. The 
scope for extra out-of-phase peaks from subdominant isocurvature modes is also 
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curtailed. The EE and TE measurements of CBI, DASI and BOOMERANG 
are mutually consistent, and, taken together rather than singly, give enhanced 
leverage for these tests. 

Subject headings: cosmology, cosmic microwave background, polarization 



Introduction 



Polarization of the cosmic microwave backgro und (CMB) at the ~ 10% level has been 
forecast for decades e.q.. lBond fc Efstathioul (119841 ). but only after a long experi mental strug- 



gle w as it detec ted and firmly est ablishe d, from measurements by the DASI (IKovac et al. 



Readhead et a 



Montrov et al. 



2006; 



2004H). CAPMAP (IBarkats et al. 



Piacentini et al.l 120061 ) exper- 



2003), WMAP fcoeut et al.ll2003h CBI 
20051 ) and Boomer ang (IMasi et al. booi 

iments. The CBI ( jPadin et al.ll2002l ; iReadhead et al.ll2004bl ) is a 13-element interferometer 



located at 5080 meters in the Chilean Andes operating in ten 1-GHz bands from 26 GHz to 
36 GHz. The CBI has been observing the polarization of the CMB since its inception in late 
1999, and since 2002 September has been operating in a polarization optimi zed configuration 



with 42 polarization-sensitive baselines. The first CBI polarization limits (ICartwright et al. 



20051 ) used only 12 polarization-sensitive baselines. As part of the polarizati on optimiza- 



tion, w e adopted the achromatic polarizers designed by Kovac and described in IKovac et al. 



( 120021 ) . The first CBI detections in the p olarization-optimized c onfiguration included data 
taken from 2002 September to 2004 May (IReadhead et al.ll2004bl ). In this paper, we present 
and analyze the implications of CBI polarization for our complete 2002 September to 2005 
April polarizatio n dataset. This represen ts a 54% increase in integration time over the results 
reported earlier (IReadhead et al.ll2004bl ). 



In § El we present an abbreviated description of the data and its analysis that leads 
to the compression of the data on to maps in £-space, a natural space for interferometry 
of the CMB, and further compression on to power spectra. A detailed description of the 
experiment, our analysis procedure and results of data quality tests will be given in Myers et 
al. (2006, in preparation). We use our improved EE and TE power spectra to test consistency 
of cosmological parameters with results forecast from TT, for minimal inflation-motivated 
tilted ACDM models in § 13 .21 and hybrid models with an additional subdominant isocurvature 
component in § 13.41 Special attention is paid to overall amplitudes and pattern-shifting 
parameters in § 13.31 and § 13.51 
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2. Processing of the CBI Polarization Data 
2.1. The Polarization Data and the CBI pipeline 



The CBI instrument is described by iPadin et al.l (120021 ) and the obs erving and data 



reducti on procedures used for CMB polarization studies are described in iReadhead et al. 



( I2004bl ). Observations were made of four mosaic fields, labeled 02 h , 08 h , 14 h , and 20 h by right 
ascension, each having roughly equal observing time. The 20 h field was a deep 4.5° x 0.75° 
20 h strip of 6 pointings separated by 45'. The fwhm of the CBI primary beam is 45'. The 
other three mosaics were 6x6 pointings each covering a 4.5° x 4.5° s quare. The CBI recorded 



visibilities on 78 baselines (antenna pairs). We note that the data in lReadhead et al.l (l2004bl ) 



used only 12 of the 13 antennas, or 66 baselines, due to a software error, which has now been 
corrected. This error resulted in a 15% loss in sensitivity, but did not bias results. The data 
were calibrated by refer ence to standard sources with an uncertainty of 1.3% in flux density 



( IReadhead et al.l l2004bl ) . Six of the antennas were set to receive right circular polarization 
(R) and seven left (L) , so each visibility measurement represents one of the four polarization 
products RR, RL, LR, LL. The copolar products RR and LL are sensitive to total intensity 
or brightness temperature T (under the assumption that the circular polarization V is zero, 
as expected), while RL and LR are sensitive to linear polarization which can be divided into 
grad-mode E and curl-mode B components. In the small angle approximation, where the 
celestial sphere can be described by a tangent plane ("image" or "sky" plane), the angular 
spherical harmonic multipoles defining the CMB radiation field, labeled by (I, m) , become 
Cartesian components of £, (£ u ,£ v ). In this same approximation, interferometer visibilities 
sample the Fourier transform of the sky brightness, with (u, v) as the conjugate variables to 
angles on the celestial plane; this u uv -space" is related to "£-space" by (£ u ,£v) = 27r(u,v). 
Each interferometer baseline is therefore sensitive to emission on angular scales centered on 
spherical harmonic multipole I = 2ttx where x is the antenna separation in wavelengths. 

The visibilities are processed by convolution with an £-space gathering kernel to pro- 
duce gridded estimators Aj<g for polarizations i =T,E,B and covariance matrix elements 
Cn{u')qq'i Ct(u')qq', Cp(u')qq' for (Gaussian) instrumental noise, (Gaussian) CMB signals, 



and projection templates ass ociated with point sources and ground spillover (IMyers et al. 



20031 : IReadhead et alJliooibh . This gridding compresses the ~ 10 7 visibilities to ~ 10 4 



£-space "pixels" (labeled by Q) without loss of essential information. 

For calculation of angular power spectra, the gridded estimators and covariance matri- 
ces are passed to a maximum likelihood procedure which estimates the CMB polarization 
bandpowers qxb in band b, X=(w')=TT,EE,TE,BB, associated noise bandpowers qNXb, a 
Fisher (or likelihood curvature) matrix Fxb,X'v whose inverse encodes the variance around 
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the maximum likelihood, and £-space window functions Wxu- The qxb are defined by the 
expansion of the spectra Cf = Y^blxbCf,- (Here Cf = £{l + l){an m a*, em ) /2n, where au m 
denotes the multipole coefficients of the signal Sj, i =T,E,B.) To determine bandpowers, 
we use (theory-blind) flat Cf shapes with top hat binning (unity inside and zero outside of 
the band). The window functions from the pipeline convert theory power spectra Cf into 
bandpowers qxb = ^^CfWxbt (^ + l/2)/(^(^ + 1)) to compare with the observed ones. 



2.2. The Power Spectra 



Fig. [TJ shows the CBI maximum likelihood bandpowers and their inverse-Fisher-matrix 
errors. The numerical values are given in Table [TJ To minimize band-to-band correlations 
only 7 bands are shown for EE, but we use about twice as many bands for our cosmic 
parameter analyses. A nonlinear transformation to a Gaussian in the "offset lognormal" 
combination \n(q X h + QNXh) is used t o give a more accu rate representation of the bandpower 
likelihood surface (IBond et al.ll2000l ; ISievers et al.ll2003l ). Fig. [2] shows the spectra of pairs of 
fields of the CBI data (which have better sensitivity than individual fields). These demon- 
strate that the remarkably good agreement between the CBI EE spectrum and the fiducial 
model for EE evident in Fig. [TJ is due to random chance. 

If the overall amplitude of Cp is set to a large value then the "nuisance" modes rep- 
resented in the construction of the matrix are projected out. T his is essential fo r groun d 
subtraction and point source projections for T. As discussed in iReadhead et al.l (j2004bl ). 
where we project out the brightest ~20% of NVSS sources in E and B, we see no evidence 
for any polarized point sources in the CBI data. In particular, for CBI the EE spectrum 
is strongly detected while the BB spectrum is consistent with zero, whereas uncorrelated 
polarized point sources give rise to (roughly) equal amounts of power in EE and BB. Neither 
polarization power spectrum rises like £ 2 as would be expected from an appreciable point 



Table 1: CBI 7-Band Power Spectra (Cg in /zK 2 ). 



£-range 


TT 


EE 


TE 


BB 


< 600 


2971 ± 260 


12.5 ±3.9 


-16.9 ±24.4 


0.8 ±2.8 


600 - 750 


1925 ± 252 


38.3 ±7.6 


-28.5 ±31.3 


-1.9 ±4.0 


750 - 900 


2475 ± 304 


3.6 ±9.9 


-77.6 ±41.9 


3.8 ± 11.0 


900 - 1050 


1126 ± 248 


47.2 ± 20.6 


-35.5 ±53.7 


-2.6 ±17.2 


1050 - 1200 


1256 ± 239 


11.3 ± 17.8 


-82.8 ±45.7 


21.1 ± 16.1 


1200 - 1500 


841 ± 137 


25.0 ± 15.2 


-61.7 ±35.7 


-9.9 ±12.0 


> 1500 


256 ± 118 


-20.3 ±35.4 


42.0 ± 59.6 


-22.6 ± 31.6 
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source signal. The spectra with and without the brightest ~20% of NVSS sources projected 
out in polarization are very similar, with no systematic trend in the differences. Finally, 
there are no sources visible in the CBI polarization maps. Consequently, in this work, we do 
not project out any sources in polarization. The nature of the 30 GHz polarized radio source 
population is currently poorly known, but some models have been devised based on extrapo- 
lation from lower frequency surveys. These generally indicate that contamination of the EE 
and BB power spectra should be negligi ble for the CBI £ range, observing frequency, and 
sensitivity levels, e.g., iTucci et al.l (120041 ) Further discussion of this issue will be presented 
in Myers et al. (2006, in preparation). 



2.3. Raw Maps and Signal Reconstructed Maps 



The gridded estimators allow us efficiently to reconstruct the polarization of the CMB. 
We use the fact that E and B are real field components that completely describe the polar- 
ization field to present E and B sky images in Figs. E](a), 0(a). These images are created by 
Fourier transforming the ground- filtered £-space estimators Aeq and Abq. These Aeq and 
Abq in turn are convolved representations of the true l-space E and B w hich are related 
to the linear polarization Stokes parameters Q and U (IKovac et al.l 120021 1 . These novel E 
and B images and £-space maps differ from the traditional polarization "headless" vector 
plots based on decomposing the images in Stokes Q and U into E-like and B-like compo- 
nents, as they are direct representations of the E and B fields and their transforms. Note 
that due to the non-local tr ansformation relat ing E and B to the Q and U fields on the 
sky, e.g-. lLewis et al.l (120021 ); iBunn et al.l (120031 ). the E and B images do not show polarized 
"objects" but coherences in the polarization field. 

The dominant ground contamination has been removed from the estimators by forming 
C 'n(C jv + C p) -1 A in the large Cp limit. Here Cjv serves to regularize Cp. Other regularizers, 
such as Cn + Ct, give virtually identical results. The images generated from the estimators 
include the effects of the observing strategy and of the gridding process, which include the 
mosaic pattern, uv coverage, the primary beam, and noise-weighting, similar to a standard 
interferometric "dirty map". Hence the raw images are not faithful reconstructions of true 
intensity and polarization. Furthermore, the Fourier transform of the regular mosaic pattern 
introduces "sidelobes" in the £-space map made from the raw estimators. Therefore, filtering 
and deconvolution can be beneficial to our maps. 

Figs. l3|5T b) show optimal Wiener-filter £-space maps. These are smoothed mean signals 
given the observations, 



;*ia> 



tot 



Cat + C p + Cj 



1 



- 6- 



Here the matrix R takes signal-space into data-space, with A = Rs + n, where s is the signal, 
B is a smoothing kernel, and n is the map-noise, including ground and source projection 
terms. One is free to choose the basis in which to describe the true sky signals s$, e.g., as 
Stokes parameters on the sky relative to a fixed sky-basis or as a set of a^ m coefficients or 
(in the small angle limit) as the set of related Fourier transform coefficients. For the CBI, 
and interferometers in general, it is natural to express the Sj in the £-basis, as in Figs. [3|5T b). 
but we also show the sky plane version in Fig. HI The = QxbCxb use d in the signal 
reconstructions are derived from the data using the measured bandpowers qxb (set to zero 
if negative) and the top-hat band matrices C^ b , and hence are "theory-blind". 

Typically the signal-space dimension would be larger than the data-space dimension so 
R would not be invertible. By modeling the (£ u , 4,)-plane only at points where we have 
estimators, the two dimensions are equal, R becomes square, and is, in principle, invertible. 
In practice R has an enormous condition number, so we remove poorly measured modes with 
eigenvalues < 10~ 4 of the largest eigenvalue of R in constructing the R _1 operator. This 
makes the reconstruction better conditioned and the results are insensitive to variations in 
the eigenvalue cut. The remaining "noise" in R _1 is controlled by reconvolution with a 
smooth regularizer B. 

We have freedom in the choice of the smoothing kernel B. The one we choose for the 
images is a natural one associated with the map. Letting the pixel be denoted by the vector 
Q = (£ u ,£ v ), the matrix R has components Rq,q+s£, where the vector 6£ goes over the 
region of £-space that contributes to the pixel in question. R tends to be only roughly 
independent of Q. We take our smoothing kernel Bsi to be the average of Rq^q+si, over all 
pixels Q. As can be seen in the insets in Figs. l3f5T b). this is very nearly the Fourier transform 
of the mosaic pattern on the sky, as is required to make the sky plane images reflect the 
area actually observed. For the strip Bse is quite asymmetric, as shown in Fig [3J For the 
other three square CBI polarization fields, the Bs£ is nearly symmetric, as shown in Fig 
There are low-level "sidelobes" in the i u and l v directions due to the mosaic spacing of 45'. 
Figure 0] shows the sky plane (as opposed to £-space) representation of the reconstructed 
signal maps for the deep strip. The agreement between the raw E image in Fig. [3](a) and 
the reconstructed E image in Fig. @] is quite good. 

To assess how well the mean field (s|A) describes the actual distribution of signal on 
the sky, it is important to see how large the fluctuations 5s = s — (s\A) are about it. For 
Gaussian signals the statistics of the (smoothed) 5s are fully described by the constrained 
correlation matrix 

(<W|A) = B{w s l - Ws^C'lRw^ (2) 
= B{w s + R\C N + Cp^R)- 1 ^ . 
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Here w^ 1 is the correlation matrix of s, the unconstrained signal in signal-space. The signal 
variance in data-space is therefore Ct = Riu^RL In terms of ws, the unsmoothed mean 
field is (s\A) = Wg R^C^A. These equations apply to a general R, e.g., when the signal- 
space and data-space dimensions are unequal. If R is an invertible square matrix, then the 
smoothed mean field is equation([T|) and smoothed realizations of the fluctuations are of form 

6s = (6s6si\A)*g, (3) 
(6s6s1\A) = B R-\C T - C T C- ]C T )YC- l B ] . (4) 

Here giQ are independent Gaussian random variables with unit variance and {5s5s^\A)z is 
the matrix square root. Equation (j3j) shows the fluctuations go to zero in modes in which 
the generalized noise is small, but approach pure signal realizations in modes in which it is 
high. Figs. EH [5](c,d) show a few examples using the Ct derived for the four CBI fields. These 
illustrate that the reconstructed signal for the deep strip is better determined than for the 
mosaic fields. 



We have also used a modified version of the CLEAN deconvolution algorithm (IHogbom 



19741 ) to do the effective inversion of R: we find the largest signal among the uv estimators 
AiQ] we place a ^-function in £-space there that zeroes out the Ajg; we subtract from each 
estimator its response to that signal; we then repeat, ending only when ~ 10 4 components 
have been found. This leads to an error in the residual less than 10~ 5 of the power in the 
original. This method has several nice features: since each estimator has compact support 
in the uv plane, the process is quite stable because the addition of a model component 
only affects a few estimators; since the model is in the same space as the data, no (time- 
consuming) Fourier transforms or (expensive) decomposition of R are needed. As in the 
standard interferometer imaging application of the CLEAN algorithm, a smooth restoring 
convolution kernel B is required to turn the set of ^-functions in £-space resulting from the 
CLEANing into a smooth transformable map, otherwise the image would have artifacts. 

We find that the results using eigenvalue cuts or this CLEAN method give very similar 
maps. The sky plane and £-space total intensity images of the 02 h CBI mosaic which we 
display in Fig [6] were constructed using this CLEAN algorithm. As expected given the 
relatively small errors on the T bandpowers, all four fields show strong T detections. 
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3. Parameterized Polarization Phenomenology 
3.1. TT, EE, TE and BB Bandpower Data 



In our parameter determinations we consider five combinations of bandpower data: (1) 
CBI TT+EE+TE bandpowers obtained from the analysis described here of the 2002-2005 
data with a bin width A£ « 75, more fine-grained than those of Table [1] (available as an 
on-line supplement); (2) CBI I = 600 t o I = 1 960 TT bands from the combined mosaic and 
deep field analysis of iReadhead et al.l f 2004a); (3) TT and TE WMAP1 bandpowers from 



the first year W MAP data 
dure described in 



jLeitch et al 



2005; 



200 



Jones et al 



(IBennett et al.ll2003l ). adopting the likelih ood mapping proce 



Verde et all j2003h : (4) DASI TT fcovac et al-lbooj ) and 3-ye ar EE+TE 



res ults; (5) the recent Boomerang B03 TT+EE +TE results (IMasi et al. 



20061 : iMontrov et al.ll2006l : IPiacentini et all 120061 ). with I < 300 TT band 



powers excluded because of overlap with WMAP (although this has no quantitative impact). 
We also omit CBI TT results for t < 360 both because of the overlap with WMAP, and 
because the (very limited) sensitivity there is coming from the sidelobes of the primary 
beam which are extremely diffi cult to measure accur ately. The Markov chain Monte Carlo 
(MCMC) package CoSMoMcfl (ILewis fc Bridldl2002l ). modified to include polarization spec- 
tra, the cross correlation between TT and EE spectra, and isocurvature modes, is used to 
calculate posterior probability distributions (including priors) for cosmic parameters. 



3.2. The Basic Flat Tilted Adiabatic ACDM Model 

The simplest inflationary paradigm is characterized by six basic parameters: uj^ = Qbh 2 , 
the physical density of baryons; uj c = Q c h 2 , the physical density of cold dark matter; 9 = 
lOO^J 1 , parameterizing the angular scale i' 1 associated with sound crossing at decoupling, 
which defines the overall position of the peak-dip pattern; ln(10 10 y4 s ), the logarithm of the 
overall scalar curvature perturbation amplitude A s , the scalar curvature power spectrum 
V s (k) evaluated at the pivot point k n = 0.05Mpc _1 ; n s , the spectral index of the scalar 
perturbations, defined by V s (k) oc &™ s_1 ; and r, the Thomson scattering depth to decoupling. 
We do not consider gravitational-wave induced components. 

Table [2] shows the broad priors we have chosen for the basic parameter ranges so that 



1 None of the conclusions drawn here are affected by using the WMAP 3-year power spectrum (released 
after submission of these results). 

' http : //cosmologist . inf o/cosmomc 
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they have little influence on our results. We also impose a weak-h prior on the Hubble 
constant H = 100/ikms -1 Mpc -1 : 0.4 < h < 1. For the flat Q t ot — 1 models considered 
here this weak-h prior has little influence on the results, although some extreme models with 
high Thomson depth are excluded. The strongest prior is the flat restriction, expected in 
most inflation models. Some parameters change significantly when the flat prior is relaxed 
faond et al.ll2003l : iReadhead eiTal1l2004ah . 



We highlight distributions for two parameters, the pattern-shifting 9/9$ and q s = 
A s e~ 2r / A S Qe~ 2ro which determines the overall Cf amplitude. We normalize relative to 
#o = 1-0442 and A s oe~ 2TO = 18.0 x 10~ 10 , the best-fit values for WMAP1+CBI+DASI+B03 
TT+EE+TE (Table [2]). The near-degeneracy between A s and r is only weakly broken at very 
low I where reionization has some influence, and at higher I through (nonlinear) secondary 
phenomena such as weak lensing or the Sunyaev-Zeldovich effect. However, the q s combina- 
tion is reasonably well determined, although there are correlations with other parameters, 
e.g.,u>b,uj c and n s , especially with polarization data only. 

The relative positions of the peaks in the TT and EE spectra are "phase-locked" by the 
physics of the acoustic oscillations at photon decoupling, with the multipole of the j th TT 
peak cx and the multipole of the j th EE peak oc (J + l/2)9~ 1 . The TE cross-spectrum 
has double the number of peaks (Fig.HJ. For WMAP1+CBI+DASI+B03 TT+EE+TE, we 
find 9/9 = 1.001 ± 0.0041 and q s = 0.996 ± 0.037 For WMAP1+CBI 9/9 n = 0.999 ± 0.005, 
to be contrasted with the 1.000 + 0.005 obtained in IReadhead et al.l (l2004bl ) for this 9 . These 
results for the mean values and standard deviations are als o very similar to those ob tained 
with other CMB data combinations, e.q jBond et al.l (120031 ); iMacTavish et al.l (120061 ) . 



For WMAP1+CBI+DASI+B03 TT+EE+TE the other basic cosmic parameters, after 
marginalization, have distributions shown in Fig. [7] with median values and la errors given in 
Table [21 The best-fit parameters for the fiducial model are also given there. Using WMAP1 
TT+TE and CBI+DASI+B03 TT gives very similar results: the inclusion of the current 
high i polarization data has little impact on parameter values for this limited set. 



For CBI+DASI+B03 EE+TE, we get 9/9 = 0.986 + 0.017 and q s = 0.86 + 0.14, in good 
agreement with the TT result. (A few other parameters are also moderately well constrained, 
but most are not, as shown in Table [2j) Our 9/9q result is not affected if we relax the flat 
prior (0.988 ± 0.018). For CBI EE 9/9 is 0.986 ± 0.031 and q s is 1.27 ± 0.48. For CBI 
EE+TE 9/9 is 0.95 ± 0.027 and for DASI+B03 EE+TE it is 1.03 ± 0.030. For CBI BB we 
obtain 9/9q = 0.95 ± 0.10. This should be interpreted as essentially the limit that we would 
get from the prior probabilities alone. This shows our 9/9$ results are data-driven rather 
than prior-driven. 
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Table 2: Cosmic parameter values for the flat tilted adiabatic ACDM model 





prior 
range 


WMAP1+CBI+DASI+B03 
TT+EE+TE (best-fit) 


WMAP1+CBI+DASI+B03 
TT+EE+TE 


CBI+DASI+B03 
EE+TE 


0/0o 


0.5 to 10 


1 


1.001 + 0.0042 


0.987 + 0.017 


Vt b h 2 


0.005 to 0.1 


0.0226 


0.0232 + 0.0013 


0.018 ± 0.005 


n c h 2 


0.01 to 0.99 


0.117 


0.114 + 0.011 


0.119 ± 0.034 


T 


0.01 to 0.8 


0.105 


0.149 + 0.086 


0.33 ± 0.18 


n s 


0.5 - 1.5 


0.960 


0.978 + 0.039 


0.92 ± 0.23 


ln[10 10 A s ] 


2.7 to 4.0 


3.09 


3.18 + 0.16 


3.37 ± 0.35 


q s = A s e- 2T /A s0 e- 2T ° 




1 


0.992 + 0.037 


0.86 ± 0.14 


n A 




0.714 


0.733 + 0.054 


0.58 ± 0.25 


Age(Gyr) 




13.6 


13.5 + 0.26 


14.4 + 0.80 






0.286 


0.267 + 0.054 


0.42 ± 0.25 


OS 




0.83 


0.848 + 0.063 


0.94 + 0.21 






12.5 


15.1 ± 5.3 


32 ± 15 


Ho 


40 to 100 


70.0 


72.6 + 5.6 


64 + 15 



The first group shows the six independent (fitted) parameters, the second group shows 
parameters derived from them. Mean values and standard deviations are given for TT+EE+TE 
data in column 4 and for EE+TE in column 5. The ranges for the uniform weak priors we 
imposed for the MCMC runs are given in column 2. The best-fit model parameters defining our 
"fiducial model" are shown in column 3. For this model 6q = 1-0437 and A S Qe~ 2T ° = 17.9 x 10~ 10 . 
Here These are sligh tly different than the parameters defining the WMAP team's best-fit 
dSpereel et alJl2003h using WMAP1 TT+TE + ACBAR TT + an earlier version of the CBI TT 

data (jPearson et alJbood ) and different priors: Q. b h 2 = 0.0224, n c h 2 = 0.111, 

n s (k = 0.05) = 0.958, r = 0.11, Hq = 72. This was the fiducial model used in lReadhead et al. 



0.05) 

(|2004bl ). Fig. [1] shows that the two are very similar visually. 
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3.3. A Peak/Dip Pattern Test 

Single-band or broad-band results using theoretically motivated Cf h shapes can also be 
produced by our pipeline. These allow complete mapping of the full likelihood surfaces 
without using the compressed bandpowers. The one-band model Cf = q s Cf , with Cf^ 
the fiducial adiabatic model, yields for CBI EE q s = 1.02±0.14 (68%) and a 11.7 a detection 
relative to q s = 0; CBI EE+TE gives a 12.4 a detection. These can be compared with the 



6.3 a DASI EE detection re ported in iLeitch et all (12004? ) and the 8.9 a CBI EE detection 



reported in lReadhead et al.l (l2004bl ) (with no polarization point sources projected out, 7.0 a 
with 20% removal). Alone, the new CBI TE data give q s = 1.02 ± 0.24 and a 4.2 a detection 
relative to q s = 0. The CBI TT data yield q s = 1.12 ± 0.05 which is a 95 a significance 
detection versus q s = 0. 

To further complement the MCMC determinations, we consider the two-parameter tem- 
plate modejfl Cf = q s CfJj^, e y evaluated on a grid in (q s , 9/9 ). The other cosmic parameters 
are fixed at the fiducial model values. We restrict 9/9 to lie between 0.3 and 1.7, the range 
of our grid. Fig. Ufa) shows how the EE peak/dip pattern shifts for the polarization. The 
9/9 -q s likelihood contours in Fig. [8^b,c,d) show that for each of the EE polarization datasets 
there is a multimodal probability structure. For example, for CBI, apart from the 9/9q « 1 
solution, there is another with the third polarization peak shifted and scaled to fit the second 
peak of the fiducial model. There is a strong probability minimum in between the two. This 
multiple solution disappears when DASI and B03 are combined with CBI, yielding the well- 
determined 9/9 = 0.988 ±0.018, q s = 0.97 ±0.09 for EE±TE. These are in good agreement 
with the 0.986±0.017 and q s = 0.86±0.14 MCMC numbers determined by marginalizing over 
the other five cosmic parameters. The multimodal aspect is strongly suppressed in MCMC 
just because of the extremely weak-h prior we impose, but correlations among parameters 
lead to the larger errors in q s . For EE alone, the template grid gives 9/9q = 0.993 ± 0.027, 
q s = 1.00 ± 0.11 and the marginalized MCMC gives 9/9 = 1.003 ± 0.027, q s = 1.07 ± 0.30. 

Fig. fS] shows the best-fit EE+TE power spectrum. It looks remarkably like the TT 
fiducial model forecast. 



5 In 



Readhead et al. ( 2004b| ) we also described a 2-parameter "sliding comb" test of the phase-relationship 



between TT and EE. This involved an underlying smooth Cf with a sinusoidal pattern characterized by 
an angular phase shift cf> designed to give the fiducial model forecast for EE when (j> = 0. The best-fit CBI 
EE phase was 21° ± 40° with amplitude q s = 1.07 ± 0.21; the new data give 13° ± 36° and 1.05 ± 0.13. 
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3.4. Adding a Subdominant Isocurvature Mode 



Isocurvature modes which could lead to a measurable signal in the CMB may arise in 
multiple scalar field models during inflation or they can be generated after inflation has 
ended. Necessary ingredients to impact CMB and large scale structure (LSS) observations 
include: association with a component of significant mass-energy, such as baryons, cold 
dark matter, or, possibly, massive neutrinos (hot dark matter); sufficiently large primordial 
fluctuations in the entropy-per-baryon, the entropy-per-CDM-particle or the entropy-per- 
neutrino. A concrete CDM realization is the axion. For cosmic-defect-induced isocurvature 
perturbations, which could also arise near the end of or after inflation, the mass-energy is in 
the defects. 

Isocurvatur e perturbations from quantum zero point fluctuation s in inflation have a well- 
defined pattern fafetathiou fc Bondlll986l : lBond fc Efstathioulll987h in TT, EE and TE, with 
no BB predicted (except through lensing). The peaks and dips are predicted to be out of 
phase with those from adiabatic modes, as shown in Fig. (For defects it is difficult to get 
any peaks and troughs at all.) Further, there is a large "isocurvature effect" predicted at 
lower I relative to that at high i where the peaks are. 

A pure isocurvature mode does not fit the current TT data unless the primordial isocur- 
vature power spectrum V- mo {k) is designed to mimic the observed Cj T pattern with its own 
peak/dip structure and overall fc-dependent blue tilt. Although highly baroque in terms of 
inflation models, such radically broken scale invariance is possible for isocurvature perturba- 
tions just as it is for the adiabatic case. Polarization (and LSS) data help by breaking such 
severe degeneracies with the cosmic parameters. 

A detailed analysis of a general set of isocurvature initial conditions for f our cosmological 
fluids using CMB and LSS data has been undertaken in lMoodley et al.l (120041 ) . If one includes 
all allowed isocurvature and adiabatic perturbations, and correlations between them, the 
current CMB and LSS data still allow a substantial amount of isocurvature perturbations. 
However simpler and more realistic models that only include an isocurvature perturbation 
in one fluid are more strongly constrained. 

Here we assume Gaussian-distributed CDM isocurvature perturbations and add two 
extra parameters beyond our base adiabatic six: two amplitude ratios, Ri = Vi SO (ki) /V s (ki) , 
at two pivot wavenumbers k iy one at small scale, k 2 = k n = 0.05Mpc _1 and one at large 
scale, k\ = 0.005Mpc _1 . A (constant) primordial spectral index defined by Vi SO (k) oc k niao 
follows: n iso = n s — 1 + \xi(R 2 / Ri) / hi{k 2 / k{) . 

We find that for neither Rj pa r amete r is there evidence for an isocurvature detection, 
in agreement with iMacTavish et al.l (120061 ) who used the same R\-R 2 parameterization. We 
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find for WMAP1+CBI+DASI+B03 TT+EE+TE, 95% confidence upper limits of R x < 0.26 
and R 2 < 1-7 on the higher wavenumber scales which CBI probes. This translates into 
steeper n iso values being more allowed than the n iso « nearly scale invariant ones. The 
CBI EE+TE data only limits R x < 18 and R 2 < 54, whereas CB1+DAS1+B03 EE+TE 
gives Ri < 9 and R 2 < 30. 

Inflation models more naturally produce nearly scale invariant isocurvature spectra, 
with n iso « 0, just as one often gets n s ~ 1 for adiabatic perturbations. The tilts from 
theory are also more likely to be red (n iso < 0) than blue (n iso > 0). However the data more 
strongly constrain red models than blue. 



3.5. Constraints on Interloper Isocurvature Peaks 

To focus attention on the high i polarization results, we now fix n iso to be the extremely 
blue 3, the white noise 'isocurvature' spectrum, with no spatial correlation. Large angular 
scales in Q^ %s °) are highly suppressed and the isocurvature peaks and troughs emerge looking 
somewhat like an ^-shifted version of the adiabatic spectrum, as shown in Fig. [9j (The 
which looks even more like a shifted version of the fiducial model, gives similar 



results to those given here. See iMacTavish et al.l (120061 ) for a treatment of both n iso = 2 and 



3 cases.) Although Vi SO {k) is so steep for such blue spectra that it must be regulated by a 
cutoff at high k >> k 2 , Cf^ ^ has a larger natural damping scale so we do not need to add 
another parameter. 

The 2-parameter template model, Cf = qgCf^ + qisoCf^ ^, therefore tests at what level 
an interloper set of isocurvature peaks would be allowed by the CMB data which, as we have 
seen in the 6/6q test, prefer the adiabatic peak positions. We normalize qi so so that qi so = 1 
corresponds to the same power in £ EE ( lso ) as j n cf E ^ over a band from I = 400 to 1200. We 
find q iso ~ R 2 q s /80. Fig. Mjo) shows a strong preference for the pure adiabatic mode and no 
isocurvature detection, with q iso = -0.01 ± 0.21 for CBI EE, 0.24 ± 0.15 for CBI EE+TE 
and 0.10 ± 0.11 for all of the polarization data. 

We also let the full 7 cosmological parameters vary, using CosmoMC to evaluate the 
probability distribution for R 2 . This is a different exercise than the 2-parameter case: to 
match the data, the other parameters are adjusted by CosmoMC to make the isocurvature 
troughs and peaks interfere with the adiabatic peaks and troughs, respectively, to mimic no 
interloping at all. For CBI EE+TE we get R 2 < 76 whereas for CBI+DASI+B03 EE+TE we 
get R 2 < 44. For WMAP1+CBI+DASI+B03 TT+EE+TE, we get R 2 < 3.0. We note that 

c TT(iso) /c TT(s) ^ 

same as for the EE ratio. Thus the upper limits correspond to 
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an allowed CMB contamination of this subdominant component of only ~ 3%. For EE+TE, 
Qiso — 0.16 ± 0.21 with an upper limit of 55%, similar to the template value. 



4. Conclusions 

In this paper we present the results from 2.5 years of dedicated polarization-optimized 
measurements with the Cosmic Background Imager. From this data, we estimate the TT, 
TE, EE and BB CMB angular power spectra. The EE power spectrum gives us a 11.7a 
detection of polarization, the strongest thus far, while TE is measured at 4.2<r versus zero. 
The BB spectrum gives a 95% confidence upper limit of 3.8 fiK 2 . 

We introduce a novel method for the reconstruction of £-space maps of E and B. Images 
of the E and B fields on the sky are formed by Fourier transform of the £-space maps; this is 
a new way of representing CMB pola rization and is compleme ntary to the standard Stokes 



Q and U images shown previously in iReadhead et al.l (j2004bl ). The E-mode detection and 



the lack of one in B is evident in both the raw maps and the reconstructed Wiener-filtered 
images of the 20 ft strip and is also evident in the square mosaic fields. We have also verified 
that signal- map fluctuations, shown in Fig. 0](c,d), about the mean signal in Fig. 0(b) do 
not obscure this clear detection: the 20 /l strip is indeed dominated by the CMB polarization 
signal. The signal maps of the total intensity, an example of which is shown in Fig [H] for the 
02 h mosaic, also show very strong detections. 

An analysis of a six- dimensional space of cosmological parameters shows that the pat- 
terns and amplitudes in the EE, TE and BB data are entirely consistent with the basic 
inflation-based model predictions from TT, a result considerably strengthened by the new 
CBI EE+TE data. The combined CBI+DASI+B03 EE+TE data further sharpens this con- 
clusion. This is particularly evident in Fig. [8] which shows that 6/6q, parameterizing the 
angular scale associated with sound crossing at decoupling and hence the peak-dip pattern, 
is pinned down to the value we obtain from TT alone. 

We finally explore a restricted physically-motivated class of models with combined, but 
uncorrelated, adiabatic and isocurvature perturbations. We find that there is effectively no 
evidence for an isocurvature mode in the data. Furthermore the data rule out a possible 
family of interloper peaks which would be out of phase with the standard flat adiabatic 
predictions. This strengthens our claim that cosmological models with an additional isocur- 
vature mode are disfavored by the current polarization data. 
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Fig. 1. — The blue crosses show power spectra Cf as a function of multipole I derived 
from the CBI 2002 September - 2005 April data, for the total intensity TT, the grad po- 
larization EE, the curl polarization BB, and cross-correlation TE. The magenta dashed 
curve is a flat power-law ACD M model which best-fits the WMAP, CBI, and ACBAR 
CMB data f Spergel et al. 2003 ). It is very nearly the black fiducial Cf which best-fits 
WMAP1+CBI+DASI+B03 TT+EE+TE, whose parameters are given in Tabled The black 
asterisks show the expected values of the bandpowers calculated from the fiducial model us- 
ing the window functions. The \ 2 values of the data relative to the fiducial model are 13.1, 
2.25, 2.90, and 8.43 for TT, EE, BB, and TE for 7 degrees of freedom. As expected, the BB 
spectrum is consistent with zero, with a single-band amplitude of 0.2 ± 1.6 £tK 2 , and a 95% 
upper limit of 3.76 /iK 2 . (The grey circles show the results of iReadhead et all (j2004bl ) for 
comparison. The shrinking in error bars is primarily due to the 54% increase in the data.) 
The CBI's \ow-£ (£ ^ 360) response is set by the details of the sidelobes of the primary beam, 
which are difficult to measure; analysis with finer bins shows that the apparent discrepancy 
in the first TT bin plotted here is confined to i < 360 where the wings of the primary beam 
are picking up the first Doppler peak. Consequently we neglect the I < 360 TT from CBI 
in further analysis. The data have been offset in i for clarity. 
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Fig. 2. — The power spectra for the 6 distinct pairs of CBI fields, each denoted by the colors 
and symbols indicated in the TT panel. For example, the blue +'s show the 02 h + 08 h fields 
and the blue x's show the 14 h + 20 h fields. The pairs of spectra plotted with the same color 
are disjoint data sets and hence have meaningful x 2 values. For EE, the x 2 values of the 
disjoint subsets are 2.98 (02 h + 08 h vs. U h + 20 h ), 8.55 (02 h + U h vs. 08^ + 20 h ), and 3.67 
(02^ + 20 h vs. 08^ + 14 h ) for 7 degrees of freedom. This confirms that the remarkably low 
X 2 (2.25 for 7 dof) of the joint EE spectrum relative to the fiducial model is due to random 
chance. The points have been spread in £ about the bin center for clarity. 
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Fig. 3. — (a) Raw images of the (approximate) E and B signals from the 20 h deep strip for 
£ < 1000. The circles show the fwhm of each of the six CBI pointings in the strip. A filter has 
been applied to remove the mean-mode (ground) signal, which would otherwise dominate 
the maps. No other filters are needed: the raw 20 h data are dominated by the CMB signal. 
For the data plotted here, simulations of pure-E and pure-B signals show that the mixing 
between E and B from sky coverage is ~5% in variance. The E map, which is predominantly 
signal, has a variance 2.8 times larger than that in the B map (which is consistent with the 
noise). Because Fourier transforms preserve power, the E/B variance ratio is the same in 
£ space, (b) The modulus of the optimal (Wiener-filtered) signal maps seen in the gridded 
(£ u , ^)-plane using our signal reconstruction process, where (£ u ,£v) = 2n(u,v)). The maps 
have been normalized to \fC~t, equivalent to plotting a^ i?Ti 's for a full-sky map. The radius 
of the magenta circles is £ = 630, the peak of the CBI's sensitivity. The EE to BB power 
ratio in the Wiener-filtered, reconstructed maps is 14.8. The upper inset in the B map is the 
average mosaic smoothing kernel B we chose to smooth the reconstructed maps. The strong 
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Fig. 4. — The E and B reconstructed signal sky-plane images of the 20 h deep strip, the 
transform of the £-space images in Fig. [3]^b). These should be compared with the raw images 
in Fig. |3](a). 
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Fig. 5. — (a) Raw maps of the (approximate) E and B signals from the 02 h mosaic field for 
£ < 1000. The two northern strips in this field have substantially less data than the four 
southernmost strips; to keep the noisy northern strips from visually dominating the map, 
they have not been included in the maps. A filter has been applied to remove the mean- 
mode (ground) signal, which would otherwise dominate the maps. The variance in the E 
map is 1.71 times larger than the variance in the B map (which is consistent with the noise). 
The circles show the fwhm of each of the 24 CBI pointings used in the maps (all 36 mosaic 
pointings are used in producing the spectrum), (b) The modulus of the optimal (Wiener- 
filtered) signal maps seen in the gridded (£ u , £ v ) plane using our signal reconstruction process, 
where (£ u ,£ v ) = 2ir(u,v)). The maps have been normalized to y/Ci, equivalent to plotting 
a^ m 's for a full-sky map. The radius of the magenta circles is £ = 630, the peak of the 
CBI's sensitivity. The striping apparent in the E image is due to the ground filter removing 
Fourier modes equal to the separation of the mosaic fields in RA. The EE to BB power ratio 
is 24. The ring at £ ~ 1000 in the B image is due to the ~ 1.3cx bandpower value in the 
fifth BB bin which enters into the filter. The upper inset in the B map is the average mosaic 
smoothing kernel B we chose to smooth the reconstructed maps. The strong anisotropy in 
£ u , £ v in Fig. [3] is not there in this square map. The compactness of B relative to that of Fig. [3] 
demonstrates the improved £-space resolution of a mosaic, at a cost of signal to noise ratio. 

TVio nrncc naf^rn ^vf^nrlincr frnm i~)ne> r'tm^ral Ti^aV nf t% aricpc frnm f Tip* Pniinpr ("rfincfrirm nf 
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Fig. 6.— (a) The reconstructed sky image for the Wiener-filtered 02 h mosaic in total 
intensity. Both ground and sources (which would otherwise dominate the map) have been 
removed as part of the filter. This image has been produced using the uf-plane CLEAN 
procedure described in § 12.31 The circles show the fwhm of each of the 24 CBI pointings 
used in Fig. [5j . (b) The modulus of the optimal (Wiener-filtered) map of (a) in the gridded 
(£ u , £ v ) plane. As in Fig. 0(b), the vertical striping is due to the ground filter, (c) A rescaling 
of (b) that brings out the structure in the £ u , 4,-plane by multiplying each pixel by I. This 
is equivalent to using £{£ + l)C^/27r instead of Ce to plot the angular power spectrum. In 
this representation, one can see hints of multiple rings of higher signal that correspond to 
the Doppler peaks in the power spectrum. 
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Fig. 7. — One-dimensional likelihoods for the cosmological parameter indicated, marginalized 
over the other parameters, (a) shows CBI EE (red), CBI+DASI+B03 EE+TE (blue). The 
black line is the prior, calculated using a parameter run with CBI BB only, (b) shows 
WMAP1 only TT+TE (black), WMAP1 TT+TE + CBI+DASI+B03 TT+EE+TE (red) 
and CBI+DASI+B03 EE+TE (blue). 
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Fig. 8. — (a) shows the polarization data from CBI, DASI, and B03, along with the fiducial 
EE prediction (black dashed), the same model shifted by 20% in 6 (black dotted), and 
the best-fit CBI+DASI+B03 EE cosmology (magenta). The polarization data pick out 
the same sound-crossing angular scale as the TT data do, with 0/9q = 0.986 ± 0.017 when 
marginalized over the other 5 cosmological parameters, (b) shows 1 and 2 a contours derived 
from CBI EE (green), DASI EE (red), B03 EE (blue) and CBI+DASI+B03 EE (magenta) 
for the 2-parameter q s and 9/9$ template model. Marginalization yields 9/9 = 0.993 ± 
0.027, q s = 1.00 ± 0.11. The yellow star marks the expected result from the fiducial model. 
The lower inset shows the same with EE+TE instead of EE only. The upper inset shows 
the sequence CBI EE (green), CBI+DASI+B03 EE (magenta), CBI+DASI+B03 EE+TE 
(brown). Marginalization yields 9/9 = 0.988 ±0.018, q s = 0.97±0.09. Grey contours denote 
CBI+B03 TT. 
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Fig. 9. — (a) c™( lso ) anc i cJ E ( lso ) power spectra for the n iso = 3 white noise isocurvature 
model (red dashed) are compared with Cf and Cj E<s ' 1 for the best-fit adiabatic fiducial 
model (black). These are the template spectra used for the 2-parameter q s -qi SO test. qi SO is 
normalized to give the fraction of the expected EE power that is observed in the data over 
the range 400 < £ < 1200. The CBI (blue asterisks), DASI (red x's) and B03 (black circles) 
EE and TE data are also shown, (b) The likelihood surface for the 2-parameter q s -qi SO model, 
for the CBI EE (light green), CBI EE+TE (dark green), CBI+DASI+B03 EE+TE (brown), 
and CBI+B03 TT (grey) data. The data strongly prefer the adiabatic over the isocurvature 
spectrum. Marginalization over the 2D distributions yields: for CBI EE q s = 1.05 ± 0.22, 
q iso = -0.01 ±0.21; for CBI TE q s = 0.81 ±0.24, q iso = 0.49 ±0.26; and for CBI+DASL+B03 
EE±TE data q s = 0.90 ± 0.10, qi SO = 0.10 ± 0.11. The polarization data are consistent with 
a single component adiabatic inflation model. 



